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More than 450 extrasolar planets are known to date. To detect these intriguing objects, many photometric 
and radial velocity (RV) surveys are in progress. We developed the Keplerian FITting (KFIT) code, to 
model published and available RV data. KFIT is based on a hybrid, quasi-global optimization technique 
relying on the Genetic Algorithms and simplex algorithm. Here, we re-analyse the RV data of evolved K3III 
star HD 240210. We confirm three equally good solutions which might be interpreted as signals of 2-planet 
systems. Remarkably, one of these best-fits describes long-term stable two-planet system, involved in the 
2:1 mean motion resonance (MMR). It may be the first instance of this strong MMR in a multi-planet 
system hosted by evolved star, as the 2:1 MMR configurations are already found around a few solar dwarfs. 

I— 8 

Introduction 

en 

,_, Most of the discovered exoplanets are found around main sequence stars, perhaps because the determination 

PLh °f stellar parameters is much easier than, for instance for giant or active stars. The formation theory of planets 

["tI hosted by such stars is still developed and not understood well. Recently, there are only known about of 

• 10 exoplanets orbiting stars with masses greater than 2 M Q (see, e.g., Extrasolar Planets Encyclopaedia, 

http://exoplanet.eu). Some recent RV surveys focus on detecting planets around the red giants, which 

I 1 are evolved main sequence stars. Unfortunately, the giants and sub-giants are difficult targets for the RV 

O technique. Usually, they are chromospherically active, pulsating, surface-polluted by large spots, and rotating 

|~) slowly. The giants produce small number of sharp spectral lines and their chromospheric activity and spots 

c/2 may change the profiles of spectral lines. This intrinsic RV variability (also known as stellar "jitter") is 

, C3, significantly larger than instrumental errors, and may be ~ 20 — 30 m/s and larger. Actually, the stellar 

activity may even mimic planetary signals (see, e.g. [1]). Hence, when interpreting the RV variability, we 

t— I may expect that different uncertainties and errors (usually, of unknown origin) may significantly shift the 

^- best fits from the "true" solutions in the parameters space. In such a case, possibly global exploration of 

Q\ the parameter space and the dynamical stability of multiple systems as an additional, implicit observable 

CN may help us to correct and verify the derived best-fit solutions for these factors, and to conclude on the 

C^ architecture of interacting systems, even when limited data are available |3]. 

O Keplerian model of the RV and optimization method 

We recall the kinematic RV model for the A^-planet system, as the first order approximation of the TV-body 
£> model: 

• rH N 

X V r {t) = y K 4 [cos{uji + !/(*)) + d COs(Wi)] + V , (1) 

^ 

where, for each planet in the system, K{ is the semi-amplitude of the signal, e« is the orbital eccentricity, Ui 
is the argument of pericenter, Ui(t) is the true anomaly, which depends on the orbital period Pi, the time of 
periastron passage Tj and eccentricity, and Vq is a constant instrumental offset. The iV-planet system is then 
characterized by N p = 5N + 1 free parameters, to be determined from observations. Let us note, that the 
inclinations and longitudes of nodes are not explicitly present in Eq. (fiT) and cannot be determined directly 
from the RV data alone, at least in terms of the kinematic model in Eq. (uB. 

To fit model Eq. (nj) to the RV data, the Gaussian Least-Squares method is commonly used. To estimate 

the best-fit model parameters, we seek for the minimum of the reduced x 2 function, minyxj(p), which 



YSC 17 Proceedings of Contributed Papers A. Rozenkiewicz, K. Gozdziewski & C. Migaszewski 

is computed on the basis of synthetic signal V r (ti,p), where ij are moments of observations V°\ s with 

uncertainties (7j (i = 1, . . . , Nrv is the number of data), and p = (K p , P p , e p , uj p , t p , Vq), and p = 1 = b,2 = 
c, . . . ,N are for the model parameters of the A-planet system. 

Clearly, even in the kinematic formulation, vxJ(p) ^ s a non-linear function. It is well known that it may 
exhibit numerous local extrema. Hence, the exploration of multi-dimensional parameter space p requires 
efficient, possibly global numerical optimization. In the past, we found that good results in this problem may 
be achieved by an application of the Genetic Algorithms (GAs) |2j, which mimic the biological evolution. 
GAs search for the best fit solutions (best adapted population members) to the model (to the environment) 
through "breading" an initial population (parameters set), under particular genetic operators (e.g., cross- 
over, mutation), and through a selection of gradually better adapted members. Although the GAs start 
with purely random initial "population", the search converges to the best fit solutions [2 J by deterministic 
way. Of course, the best fit solutions may be not unique. This is very common in the case of modeling 
RV data. Overall, the GAs are robust and quasi-global optimization technique, although they not provide 
efficiency and accuracy of fast local methods, like the Levenberg-Marquardt algorithm. The GAs are used 
in our hybrid KFIT code, developed for a few years now [3], which also makes use of the local, and accurate 
simplex algorithm [BJ. It helps to quickly refine the final population of the best fits, "grown" by the GAs. 

The results for the RV of HD 240210 

In a recent work, Niedzielski et al. [5] detected a planet hosted by evolved dwarf HD 240210. They found an 
excess in the rms of Keplerian 1-planet model, and interpret this as possible signal of an additional planet. 
We reproduced their 1-planet solution, see Fig. [I] {the top-left panel). The rms has significant scatter of 
39 m/s, compared to the mean instrumental accuracy of ~ 8 m/s. This rms excess is difficult to explain by 
the internal variability of the star, hence we tried to find a better 2-planet model. At first, we reproduced 
(see Fig. Fy the top-right panel) a tentative 2-planet solution illustrated in the discovery paper. However, 
we noticeihat the orbital periods are very similar ~ 440 and ~ 530 days, respectively, indicating strong 
4:3 MMR or 1:1 MMR. In such a case, the kinematic model is not proper anymore, due to significant mutual 
planetary interactions. Moreover, such configurations could be hardly explained by the planetary formation 
theory. Still, having in mind that this solution is not unique, we performed an extensive search for the 

local minima of yxj with the KFIT code. The statistics of gathered fits is shown in Fig. Bl through their 
projection onto orbital periods plane. In the range of orbital periods 6 [60, 3600] days, we found three, equally 
good best-fit models, which correspond to different orbital configurations, and may be resolved at the 2a 
confidence level. In the parameter maps, the mentioned 2-planet model is labeled as Fit I. Two additional 
fits with orbital periods ratio close to 3:2 and 2:1 are labeled as Fit II and Fit III, respectively (see Table 1). 
Their synthetic RV curves, with the RV of 1-planet model (thin curve) and observations overplotted, are 
shown in subsequent panels of Fig. [TJ Note, that the alternative 2-planet fits reduce the rms significantly, 
to ~ 25 m/s (i.e., by 1/3), consistent with a conclusion in [5]. As may be seen in Table 1, Fits II and III 
have large eccentricities. A question remains, whether inferred orbital configurations are dynamically stable. 
In fact, all these Fits I — III, transformed to osculating elements at the epoch of the first observation, lead to 
self-disruption of 2-planet systems. Nevertheless, remembering that stable solutions may still be found in 
their neighborhood, we did dynamical analysis with the ./V-body, self-consistent GAMP code |41 (which also 
relies on the GAs), trying to refine Fits I— III with the requirement of the long-term stability (the edge-on, 
coplanar models are tested). Certainly, at most one of these models might correspond to the real system. 
We did not found any stable orbits in the vicinity of Fit I. In a tiny neighborhood of Fit II, there is a stable 

configuration with yxj ~ 1-50 and an rms ~ 27 m/s which, as the direct numerical integrations show, is 
stable over 1 Gyr. The best result is found for Fit III as a stable solution, corresponding to the 2:1 MMR 

with -\/Xu ~ 1.36 and an rms of ~ 25 m/s, i.e., the same as in the kinematic Fit III. The dynamical map 
[1] around this solution (Fig. K3J the right-hand panel) reveals extended island of stability (~ 0.15 au). This 
fit has moderate semi-amplitude librations (~ 15 — 30°) of the critical angles 9\ = 2\b — X c — vd^ (around 
0°), 82 = 2Afc — A c — w c (around 180°), and #3 = w^, — w c (around 180°). The numerical integrations 
confirmed that its stability is preserved at least over 1 Gyr. The 2:1 MMR Fit III seems the most promising 
planetary model explaining the RV variability of the HD 240210. The 2:1 MMR is quite frequent in the 
sample of ~ 40 known extrasolar systems with jovian planets, because 5-6 configurations were reported (see, 
http://exoplanet.eu). Hence, this new system, which could be the first one around evolved star, is likely. 
We stress that solution III is found in relatively extended stability zone, unlike Fit II, which lies in a tiny, 



isolated area (~ 0.01 au, Figj3 
other, relatively extended stab 



the left-hand panel). These two maps almost overlap in the a c -range, hence 
e islands are rather excluded in this region. 
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Figure 1: The RV observations of HD 240210 [5] and synthetic curves for the best fit solutions found in this paper. 
The left-upper panel: the best 1-planet fit. Subsequent figures are for the best fits labeled with I. II. Ill in Fig. [2] For 
a reference, all these plots are accompanied by the model curve of the 1-planet fit. See Table 1 for orbital elements. 
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Figure 2: Statistics of the Keplerian best fits gathered with the KFIT code, projected onto the (Pb, P c )~plane of 
orbital periods. Solutions marked with white circles and labeled by I, II, III are for the best fits in different "islands" 
of the parameters space. Black filled circles are for all fits within formal 3er-level of Fit I; for a reference contours 
with 3,2,1 c-level of this solution, found through extensive, systematic scanning |4J of the (Pf,, P c )-plane with the 
Levenberg-Marquardt algorithm, are plotted as curves of increasing thickness. See Table 1 for model parameters. 
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Figure 3: Dynamical maps in terms of the MEGNO indicator [I] around coplanar, edge-on, GAMP iV-body fits (crossed 
circle): Fit II (3:2 MMR, ^xl = 1 50; rms ^ 27 m/s, the left-hand panel), and Fit III (2:1 MMR, ^xT = 1.36, mis 
~ 25 m/s, i/ie right-hand panel). White color is for unstable configurations, black is for stable solutions. Orbits at 
the epoch of the first RV of Fit III in terms of (m [nij], a [au], e, w [deg], A4 [deg]) are the following: (4.11, 1.14, 0.284, 
304.8, 101.8) b , (2.27, 1.83, 0.592, 162.3, 304.2) c for planets b, and c, respectively, V = 21.57 m/s, stellar mass is 
0.82 M 0; Fit II is: (4.50, 1.143, 0.217, 249.24, 163.87) h , (1.94, 1.506, 0.562, 53.09, 243.39) c , V = 8.98 m/s. 



Table 1: Keplerian model parameters of 2-planet best fit solutions to the RV of HD 240210. Formal measurement 

27. 



e rescaled by adding the stellar jitter of Oj = 


= 20 m/s in quadrature. To = 


= 53, 000 days, N v = N RV ~ I 


the best fit 


I 


II 


III 


parameters 


b c 


b c 


b c 


P [day] 


540 ± 29 441 ± 27 


484 ± 14 667 ± 54 


485 ± 18 994 ± 97 


K [m/s] 


131 ±29 85 ± 30 


147 ±16 82 ± 36 


129 ±30 63 ± 23 


e 


0.05 ± 0.09 0.20 ± 0.14 


0.18 ± 0.14 0.74 ± 0.38 


0.30 ± 0.13 0.58 ± 0.35 


w [deg] 


204 ± 51 356 ± 73 


261 ±41 36 ± 76 


302 ± 35 163 ± 87 


t - T [days] 


352 ± 65 610 ± 109 


484 ± 45 502 ± 72 


536 ± 31 326 ± 156 


Vb[m/s] 


-3±7 


12 ± 9 


22± 12 


/ 9 

V ^v 


1.35 


1.36 


1.36 


rms [m/s] 


25.2 


25.2 


25.2 



Conclusions 

Extrasolar planets hosted by giant or evolved stars bring important border conditions for the planet formation 
theory. In this work, we re-analysed the literature RV data for evolved dwarf HD 240210, with our KFIT 
code relying on quasi-global GAs. In the reasonable range of orbital periods less than 3600 days, we found 

three Keplerian solutions, which have the same yxj and an rms. By further dynamical analysis of these 
best-fit models, we selected the most likely, stable solution, which corresponds to 2:1 MMR, and is located in 
relatively extend zone of dynamical stability. Overall, if the 2-planet configuration is assumed, the dynamical 
constraints seem rule out other two models, but only new observations may confirm the 2:1 MMR hypothesis. 

Acknowledgements. This work is supported by Polish Ministry of Science, Grant 92/N-ASTROSIM/2008/0. 



References 

[1] Berdyugina, S. V., LRSP, 2(8) (1995) 

[2] Charbonneau, P., ApJS, 101, 309 (1995) 

[3] Gozdziewski K., Migaszewski, C, A&A, 449, pp.1219-1232 (2006) 

[4] Gozdziewski K., Migaszewski & C, Musieliflski, A., Proc. IAU Symp. 249, pp.1219-1232 (2008) 

[5] Niedzielski A., Nowak G., Adamow M. & Wolszczan A., ApJ, 707, pp. 768-777 (2009) 

[6] Nelder J. A., Mead R., Comp. J., 7, pp. 308-313 (1965) 



